Scale-free network on a vertical plane 
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A scale- free network is grown in the Euclidean space with a global directional bias. On a vertical 
plane, nodes are introduced at unit rate at randomly selected points and a node is allowed to 
be connected only to the subset of nodes which are below it using the attachment probability: 
■Ki(t) ~ ki(i)l a . Our numerical results indicate that the directed scale-free network for a — 
belongs to a different universality class compared to the isotropic scale-free network. For a < a c 
the degree distribution is stretched exponential in general which takes a pure exponential form in 
the limit of a — > — oo. The link length distribution is calculated analytically for all values of a. 
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It has been seen in many branches of Statistical Physics 
that a global directional bias in space has strong effect 
on the critical behaviour of simple models. Introduction 
of a preferred direction in the system reduces the degrees 
of freedom of the constituting elements of the system, 
which shrinks the configuration space available to the 
system compared to the undirected system. As a result 
a directed system is simpler and quite often tractable 
analytically. Examples include Directed percolation [lJL 
Directed Sandpile Model Q, Directed River networks Q 
and Directed Self-avoiding walks Q etc. 

Over the last few years it is becoming increasingly ev- 
ident that highly complex structures of many social [l| , 
biological HI]} or electronic communication H El net- 
works etc. cannot be modeled by simple random graphs. 
For example in the well known random graphs by Erdos 
and Renyi, the degree distribution P(k) is Poissonian 
(de gree A: of a vertex is the number of edges attached to 
it) [10(. In contrast, it has been observed recently that 
the nodal degree distributions of many networks, e.g., 
World-wide web and the Internet have power law 
tails: P(k) ~ k 1 . Due to the absence of a characteristic 
value for the degrees these networks are called 'scale-free 
networks' (SFN) [H El El • Barabasi and Albert 
(BA) generated scale-free graphs where a fixed number 
of vertices are added at each time and are linked with 
a linear attachment probability [llj. On the other hand 
some of these networks are directed networks whose links 
are meaningful only when there is a connection from one 
end to the other but not the opposite, e.g., the World- 
wide web H, the phone-call graph EH and the citation 
graph EH- 

However, there are networks in which the nodes are 
geographically located in different positions on a two- 
dimensional Euclidean space e.g., electrical networks, In- 
ternet or even in postal and transport networks etc. The 
edges of the graphs representing these networks carry 
non-uniform weights which in most cases are either equal 
or proportional to the Euclidean lengths of the links. In 



these networks a relevant question is how to optimize the 
total cost of the connections e.g., electrical wires, Ether- 
net cables or say travel distances of postal carriers E3- 
On the other hand a detailed knowledge of link length 
distribution is also important in the study of Internet's 
topological structure for designing efficient routing proto- 
cols and modeling Internet traffic. For example, Waxman 
model describes the Internet with exponentially decaying 
link length distribution jig. Yook et. al. observed that 
nodes of the router level network maps of North Amer- 
ica are distributed on a fractal set and the link length 
distribution is inversely proportional to the link lengths 
E3- Other models of networks on Euclidean space are 
also studied in the literature [jfl EH, E3 | . 

In this paper we studied the effect of a global direc- 
tional preference on the statistics of scale-free networks 
embedded in the Euclidean space. A typical link in this 
model must have a positive component along some pre- 
ferred direction. Similar to the directed versions of well 
known models of Statistical Physics P, II II 01 our spa- 
tially directed networks have different universal critical 
behaviour compared to their undirected counterparts. 

A two-dimensional network is grown whose nodes are 
the points at randomly selected positions within an unit 
square on the vertical x—y plane. To construct a network 
of N + l nodes, let (x ,xi,x 2 , -..,Xn) and (yi,U2, ■■■,Un) 
be the 2N + 1 independent random variables identi- 
cally and uniformly distributed within the interval {0, 1}. 
Let a specific set of values of the random variables 
{(x ,0),(xi,yi),(x2,y2),—,(xN,VN)} represent the co- 
ordinates of the of N + 1 randomly distributed points. 
The growth of the network starts with only one node 
(xo,0) on the bottom side of the unit square and then 
the other nodes are added one by one at unit rate ac- 
cording to their serial numbers % = 1 to N. 

We assume that the global directional bias is the grav- 
ity and acts along the —y direction which restricts the 
choice of the link: a new node can only be connected 
to a node positioned below this node. In practice when 
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FIG. 1: Pictures of the networks generated from the same 
distribution of 513 points within the unit box. A large degree 
node is visible for a DSFN in (a) and long length early links 
are observed for a DMGN in (b). 

the t-th node is introduced, we consider the subset St of 
the nodes situated below the t-th node. The t-th node 
is then connected to any node of this subset using some 
specific attachment rule. In addition, we assume that 
from each node only one link comes out but any number 
of links can terminate on this node. This condition en- 
sures that the network is a singly connected tree graph. 
Initially the 0-th point is assigned the degree fco(0) = 1. 
Link lengths are measured using the periodic boundary 
condition imposed only along the x direction because of 
the anisotropy. Depending on how a node from the sub- 
set St is selected for connection we consider the following 
two models: 

(a) Directed scale- free Network (DSFN): The t-th node 
is randomly connected to a node i of the subset St using 
an attachment probability which is linearly proportional 
to the degree kt(t) of the node i at time t as: 7ii(£) ~ fcj(t). 

(b) Directed minimal growing network (DMGN): The 
t-th node is connected with probability one to the near- 
est node in the subset St- Pictures of typical network 
configurations are shown in Fig. 1. 

A continuous tuning between these two different mod- 
els is possible by the choice of a suitable tunable param- 
eter a. This is achieved by modulating the attachment 
probability in DSFN by a link length £ dependent factor 
like: 

7T,(t) ~ h(t)£ a . (1) 

This introduces a competition between the roles played 
by the degree as well as the link length on the attach- 
ment probability. The limiting extreme cases are the 
above two models. In the case with a=0 the link lengths 
do not play any role and therefore the model corresponds 
to DSFN. On the other hand when a = — oo only the 
shortest link is selected with probability one irrespective 
of the degree of the node and therefore the model cor- 
responds to DMGN. First we study these two limiting 
cases. 




io" s 10" 4 10"' io" 2 10"' io' 
k/N 



FIG. 2: The scaled degree distribution for DSFN for network 
sizes N = 2 12 ,2 14 and 2 16 . The collapse of the data at large 
k values imply that the degree distribution exponent 7 w 2.4. 

For a scale-free network the nodal degree distribution 
has a power law tail: P(k) ~ k 1 and it obeys a finite 
size scaling form: 

P{k,N) ~ N~ v J r (k/N^). (2) 

We numerically find that the degree distribution of DSFN 
(excluding the node on the bottom line) indeed follows 
such a scaling form with rj 2.4 and £ as 1 (Fig. 2). 
This gives jdfsn — rj/C ~ 2.4. This value of TD.qfw is 
compared with 7 = 3 for the BA model of SFN [llj and 
therefore it seems that DSFN belongs to a new univer- 
sality class different from BA SFN. On the other hand 
the degree distribution for the DMGN is found to decay 
exponentially as: P(k) ~ exp(-Kk) with k « 0.74. 

For a tree graph, the branch size distribution is very 
important and the associated exponent may be used to 
characterize the graph. On a tree structure, each edge 
connects two branches of the tree. If an edge is selected 
randomly, the probability Prob(s) that any one of the 
two branches supports s nodes also decays with a power 
law tail, Prob(s) ~ s~ T and follows a scaling form: 

Prob(s) - N-^Gis/N^). (3) 

For DSFN, we obtain rjb « 2.15 and 0> ~ 1 which implies 
that tdsfn ~ 2.15 compared to its exact value 2 for 
the BA scale- free network j^. On the other hand for 
DMGN we find r\h « 2 and Cfc ~ 1 so that tdmgn ~ 2. 

The probability density distribution D(£) gives the 
probability D(£)d£ that an arbitrarily selected link has 
a length between £ and £ + d£. For the undirected scale- 
free Euclidean networks we saw that D{£) has a power 
law variation D{£) ~ £ s [HE!]. D(£) can be calculated 
exactly for both DSFN as well as DMGN in the following 
way. Let us try to assign a link to the [n + l)-th node 
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FIG. 3: Variations of the average clustering co-efficients for 
DSFN: (a) Over the whole network, C(N) ~ N~ 0N and 
(b) Over the subset of nodes having degree k only C(k) ~ 
[k(\nk)]~ Pk . Our estimates are f3 N « 0.23 and (3 k w 0.64. 



and denote y — y n +i. Let ri\ points be positioned below 
the y level and n 2 = n — m points be above this level. 

We first calculate Ddmgn(^)- The probability that 
out of ni nodes the node which is nearest to the (n + 1)- 
th node is positioned at a distance between I and I + d£ 
has two distinct contributions: One from the case of all 
y > I and the other for all y < t (due to the presence of 
the boundary at y = 0). For the first case, the probability 
of the (n + l)-th point being at a particular height y is 
given by 

ir£d£K 1=0 { n C ni }ni{y - ^/2\ n ^{\ - y) n > 
= n£d£(n)[l-Tre 2 /2] n - 1 . 

Here the weight factor from each partitioning of n into 
m and 7i2 has been taken into account. The above prob- 
ability after integration over y from £ to 1 and summed 
over all n from to oo gives the net contribution to D{£) 
for all y > I: 



n-l 



D DMG N(y>£) = n£[l - £}V%L Q n[l - ir£ 2 /2] 

Similarly the probability that the (n + l)-th point is at 
a specific height y < I is 

2f S m- 1 |rf^ i=0 {"C„ 1 }(n 1 ) 
x[y- (^sm" 1 ! +2/ ^2^2)]"i-i(i _ 
= 2£sm- 1 ^n[l - (£ 2 si n - lV j + y^J P - y 2 )]^ 1 dl. 
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FIG. 4: Variation of the exponent (5{a) characterizing the 
stretched exponential degree distribution. The continuous 
curve is a fit to the data points (circles) obtained by simula- 
tion to a form: /3(a) = a(—ct v ) — b such that (5 is extrapolated 
to zero at a c ~ —0.85. 



As before a similar sum over n and integration over y 
from y = to y = £ in the above expression gives the 
following contribution to D(£) for all y < £: 

A 



Dbmgn {y < £) 

f 1 2sin~ 1 z 

where A = I -== 

Jo sin~ l z + zyl — 



-dz. 



Hence the total distribution is given by: 

D DM gn{£) = ^{1-{1-\a)£). 
For the DSFN also one has: 

f 1 ir£ 

D D SFN(y > £) = C —dy 
Ji V 

and D DSFN{ y<£) = C f ^M** 



(4) 
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where C is a constant. The total distribution is therefore 
given by 



D DSFN (£) - C(B£-£ln(£)) 



(5) 



where B = 2 



f 1 sin 1 i 
Jo 



dz. 



From our numerical calculations we estimate C w 1.59 
and B w 1.86. For any non-zero a, the corresponding 
distribution can be obtained by simply multiplying the 
above expression by £ a . 

The clustering co-efficient C(N) of a network of TV 
nodes measures the the local correlations among the links 
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of the network. More precisely it measures the probabil- 
ity that two neighbours of an arbitrary node are also 
neighbours. If the z-th node has the degree fcj and there 
are links among the ki neighbours of i then the clus- 
tering co-efficient of the site i is: Cj = 2ej/[fci(fcj — 1)] 
whereas the clustering co-efficient of the whole network 
is: C(N) = (Cj). For a number of networks it has been 
observed that the clustering co-efficient decreases with N 
like C(N) ~ N- 0N as the network size N increases. Also 
one can define a clustering co-efficient C(k) averaged over 
the subset of nodes of degree k on the network. It has 
been also observed that C(k) ~ k^ /3k for some networks. 
We estimated these exponents for DSFN and found that 
/3jv ~ 0.23 whereas C(k) has a logarithmic modulation 
like CO) - [k(lnk)]-P« with fa « 0.64 (Fig. 3). 

Finally we study the variation of the degree distribu- 
tion P(k) with the parameter a. For finite negative val- 
ues of a the distribution fits very well to a stretched 
exponential form: P(k) ~ exp(— ck 13 ^) where (3{a) is 
expected to reach to one as a — * — oo and to zero as 
a — > a c . Our numerical estimates for (3 have been plot- 
ted in Fig. 4 with a and this data fits very well to form 
(3(a) = o(— a) u — b where the constants are estimated to 
be a w 0.47, v w 0.51 and b w 0.43. This implies that 
the stretched exponential form continues to be valid till 
a = a c where (3 = and beyond that the degree distribu- 
tion is a power law. Though from the values of a, v and 
b, a c is estimated to be —0.85 we believe a c = —1 is more 
plausible. Also our numerical results indicate that for all 
a > a c the degree distribution exponent 7 maintains its 
value of a = 0. 

To summarize, we studied the directed version of 
the Barabasi-Albert scale-free network grown on a two- 
dimensional vertical plane. Our numerical results on de- 
gree as well as branch size distributions indicate that this 
network belongs to a new universality class compared to 
its undirected version. A competition between the degree 
k of the nodes and a link length dependent factor £ a in 
the attachment probability is seen to control the network 
behaviour. In the limit a — > —00 one gets the directed 
minimally growing network with exponentially decaying 
degree distribution. However for finite negative values of 
a stretched exponential distributions are observed. The 
link length distribution is calculated analytically for all 
values of a. 
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